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o : 

In one-dimensional quantum mechanics, or the Sturm-Liouville theory, Crum's theo- 
^ ' rem describes the relationship between the original and the associated Hamiltonian systems, 

, which are iso-spectral except for the lowest energy state. Its counterpart in 'discrete' quan- 

X^/^ ' turn mechanics is formulated algebraically, elucidating the basic structure of the discrete 

quantum mechanics, whose Schrodinger equation is a difference equation. 

i | §1. Introduction 

43 : 

In the seminal paper of 1955, Crum showed, 1 ) if rephrased in the language of 
■ quantum mechanics, the existence of an associated Hamiltonian system for any given 

one-dimensional quantum mechanical Hamiltonian system under mild assumptions. 
The method or the technique is quite universal and is known under many different 
names; the Darboux transformation, 2 ) the factorisation method 3 ) or the supersym- 
CN) ■ metric quantum mechanics. 4 ) Crum himself presented his results in the traditional 

language of Sturm-Liouville systems. Since Crum's theorem elucidates the generic 
structure, many exactly and quasi-exactly solvable quantum mechanical examples 
were constructed by employing Crum's theorem and its modifications combined with 
shape invariance. 5 ) In particular, Adler' modification of Crum's theorem 6 ) is quite 
CN ' general and useful. It allows to construct an infinitely many exactly solvable poten- 

tial from any exactly solvable one. 

Recently 'discrete' quantum mechanics was introduced by the present authors. 7 )~ 9 ) 
It is a generalisation of quantum mechanics, in which the Schrodinger equation is 
a difference equation instead of differential in the ordinary quantum mechanics. In 
other words, the differential operator in the Hamiltonian is replaced by finite dif- 
ference operators, either in the pure imaginary or the real direction. This is why 
they are called discrete quantum mechanics. Many explicit examples of exactly 10 )' n ) 
and quasi-exactly solvable 12 )' 13 ) systems have been constructed. These exactly solv- 
able systems have salient features; they are solvable both in the Schrodinger and 
the Heisenberg pictures. The main part of the exact eigenfunctions are the known 
hypergeometric orthogonal polynomials of the Askey scheme. 14 )~ 16 ) On the other 
hand, the exact Heisenberg operator solutions define the creation and annihilation 
operators 10 )' n ) which, together with the Hamiltonian, form the dynamical symmetry 
algebra of the exactly solvable systems. The so-called g-oscillator algebra 17 ) is the 
most typical dynamical symmetry algebra realised in this way. 18 ) In discrete quan- 
tum mechanics, as will be shown shortly, the counterpart of Crum's theorem holds 
and it also elucidates the generic structure of one dimensional systems. As in the 
ordinary quantum mechanics, Crum's theorem and its modifications are useful for 
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constructing various exactly and quasi-exactly solvable discrete quantum mechanical 
examples. In particular, its modification a la Adler 19 ) will also provide an infinite 
number of exactly solvable ones based on any exactly solvable one. When applied to 
discrete quantum mechanics with real shifts, the modified Crum's theorem generates 
an infinite or finite series of exactly solvable birth and death processes based on 
any known exactly solvable one. 20 -* The insight obtained Crum's theorems and their 
modification, in the ordinary and discrete quantum mechanics, is essential for the 
recent derivation of infinite numbers of shape invariant systems and new exceptional 
orthogonal polynomials. 21 )' 22 ) 

In this paper, we present the discrete quantum mechanics version of Crum's 
theorem. Under mild assumptions, we construct algebraically an associated Hamil- 
tonian system for any given one-dimensional discrete quantum Hamiltonian system. 
That is, if a Hamiltonian together with the full discrete energy spectra and the cor- 
responding eigenfunctions are given, then the associated Hamiltonian together with 
the full discrete energy spectra and the corresponding eigenfunctions are constructed 
as in the original Crum's paper. 1 ) This process goes on indefinitely, since the first 
associated Hamiltonian system generates the second associated Hamiltonian system 
and so on. Due to the essential distinction between the differential and difference 
equations, several technical assumptions are necessary for the discrete Crum's the- 
orem. For example, in one-dimensional ordinary quantum mechanics, the energy 
spectra are non-degenerate and the oscillation theorem holds. That is, the n-th 
excited state wavefunction has n zeros. These two properties are not necessarily 
shared by the generic discrete QM. The hermiticity is almost trivial in ordinary QM, 
but it can only be proven after the explicit form of the groundstate wavefunction 
is obtained in discrete QM. 9 )' 12 ) Another distinctive feature is the wavefunction. In 
ordinary QM, the wavefunction is a complex valued real function, defined on the real 
line, or half line or on a line segment. In one dimension, the wavefunction can be 
chosen real. In discrete QM with pure imaginary shifts, the wavefunction undergoes 
shift operations tp(x) — ► ip(x ± ry), 7 G R. Thus we require the analyticity of the 
wavefunction with its domain including the real axis or a part of it in which the 
dynamical variable x is defined. In spite of these differences, the generation of the 
associated Hamiltonian system goes almost parallel in the ordinary and discrete QM, 
since the main part is algebraic. 

This paper is organised as follows. In section two, Crum's theorem is recapitu- 
lated in some detail in the language and notation of ordinary quantum mechanics. 
This explains the underlying logical structure of the associated Hamiltonian system, 
which is also shared by the discrete quantum mechanics version. The discrete version 
of Crum's theorem is stated and proved in section three. The final section is for a 
summary and comments. 
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§2. Ordinary Quantum Mechanics 

Let us start with a generic one-dimensional quantum mechanical system having 
discrete semi-infinite energy levels only: 

Q = S <£i <£ 2 < ■■■ ■ (2-1) 

Here we have chosen the constant part of the Hamiltonian so that the groundstate 
energy is zero. Then the Hamiltonian is positive semi-definite and in one dimension 
all the energy levels are non-degenerate. It is well known in linear algebra that any 
positive semi-definite hermitian matrix can be factorised as a product of a certain 
matrix, say A, and its hermitian conjugate A* ■ Similar factorisation applies to the 
present quantum mechanical Hamiltonian 7i, which takes a very simple form: 

w = p2 + [,(,., =p , + (_!!) + _J_1 > p = _,_ {2 . 2) 

dx dx dx dx 

Here a real function W(x) € C°° is called a pre-potential and it parametrises the 
groundstate wavefunction (po(x), which has no node and can be chosen real and 
positive: 

4> (x) = e w ^\ (2-4) 

It is trivial to verify 

AM x ) = o => nM x ) = o- (2-5) 

All the eigenfunctions are square-integrable and orthogonal with each other and form 
a complete basis of the Hilbert space: 



/ 



H<j) n {x) = SnMx), n = 0, 1, 2, . . . , (2-6) 
t> n (x)*(f>m(x)dx = h n 8 nm , 0<h n <oo, n,m = 0,1,2, (2-7) 



It is well-known that the ra-th excited wavefunction (j) n (x) has n zeros in the interior. 
For simplicity we choose all the eigenfunctions to be real. Here are a few examples: 

W{x) = -\x 2 , U(x) = x 2 - 1, -oo < x < oo, (2-8) 

W{x) = -\x 2 + glogx, U(x) = x 2 + 9 ^ 9 ~ ^ - 1 - 2g, g > 1, < x < oo, (2-9) 

W(x) = ^logsinx, U(x) = 9 ^ 9 ~ ^ - g 2 , g > 1, < x < vr. (2-10) 

sin x 

They all lead to well-known exactly solvable quantum mechanics whose eigenfunc- 
tions consist of the classical orthogonal polynomials, the Hermite, Laguerre and Ja- 
cobi polynomials, respectively. In each case, the squared groundstate wavefunction 
4>o(x) 2 gives the orthogonality weight function for the polynomials. 
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Next let us define an associated Hamiltonian Ti.^ by simply changing the order 
of A and A^. 

H^=AA^. (2-11) 
For later convenience, let us attach the superscript I ' to all the quantities in the 
original Hamiltonian system, TL^ = TL, 4>n\x) = f 4n(x), A^ = f A, U^(x) = f U(x), 
W[°](x) = f W(x). By construction A and A^ intertwine Tt^ and Tt^: 

A [o] H [o] =n [i] A io] i A m n m = h[°UW. (2-i2) 

We will show that the associated Hamiltonian system TC^ is iso-spectral to the 
original Hamiltonian system and the eigenfunctions are in one to one correspon- 
dence, except for the groundstate. Thanks to the first of the above relation (|2- 12|) . it 
is trivial to verify that the eigenfunctions of the associated Hamiltonian system Ti^ 
are generated algebraically by multiplying A^ to the eigenfunction of the original 
system: 

(a-) Wa®<$ J 4 11 = SnKSnm, n, m = 1, 2, . . . , (2-13) 

wW^(x) =£: n 4 1 ](x), ti = 1,2,.... (2-14) 

Suppose the associated Hamiltonian Ti\ 1 ^ has an eigenfunction (p'(x) with the eigen- 
value £' other than those listed above: 

H [1] 4\x)=£'4\x). (2-15) 

Again, thanks to the second of the relation (|2-12j) . it is trivial to verify 

H l0] A [0] W(x) = £'A [0] W(x). (2-16) 

Due to the completeness of the spectrum of the original Hamiltonian 7i^°\ the pro- 
visional eigenvalue £' must belong to the above spectrum (|2-14j) for n = 1, 2 In 

other words, £' cannot be vanishing, £' ^ 0. Suppose that is the case {£' = 0), then 
4' is annihilated by .A^t. Surely there exists a solution of a first order differential 
equation A^<p'( x ) = 0, 4'i x ) = 1/4$ {x) = e~ w ^ . But it is obviously non square- 
integrable and it does not belong to the Hilbert space of the associated Hamiltonian 
system. After Crum, 1 ) we can show in the following way that 4n\x) has exactly 
n — 1 zeros. Since 4$ has exactly n zeros, the relation 

4 [ n\x) = d , <$ ] (x) s 

4 [ q\x) dx ^4$\x)^ 

tells through Rolle's theorem that 4>n has at least n — 1 zeros. From the relation 

^(4 [ o\x)4^ ] (x)) = -£ n 4t\x)4^(x) 1 

we find that 4n^ has at most n—1 zeros. Thus we have established that the associated 
Hamiltonian system Ti^ is iso-spectral to the original Hamiltonian system T~6^ and 
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the eigenfunctions are in one to one correspondence, except for the groundstate 
with the wavefunction </>q°'(x). If the groundstate energy £\ is subtracted from the 
associated Hamiltonian H^, it is again positive semi-definite and can be factorised 
as above: 



H [1] = A m A [1] +£i= P 2 + U [1] (x)+£ 1 , (2-17) 

(2-18) 



dx dx dx dx 



A ll] $ ] (x) = 0. (2-20) 

As shown above, the groundstate wavefunction (x) has no node. Note that (|2- 1 1|) . 
(f2~T7j) and (12091) imply the Riccati equation for dW [1] (x)/dx: 



/dW ll] (x)y d 2 W^(x) _ /dW(x)y d 2 W(a 
\ dx J dx 2 V dx J dx 2 



Si. (2-21) 



Then by reversing the order of A^ and A^ the second associated Hamiltonian 
system can be defined. This process can go indefinitely. 

Here we list the definition of the s-th quantities step by step for s > 1, 

H [ s ]drf^ [s _ 1U[s _ 1]t + ^_ i) (2 22) 

^(x)^A^t 1] (x), (n>s), (2-23) 

e wW(x)def (2 24) 

A [s]W d <M l ' ] W A m d dWM(x) 
dx dx ' dx dx ' 

\ dx J dx z 
Then we can show the following for n > s > 0, 

wH0H(x)=^4 s ](^), (2-27) 

4> [ n ] (x) : real function, (2-28) 

A [s] ^ ] (x) = 0, (2-29) 

W W = t^M + s s = p 2 + £/M ( x ) + £ s . (2-30) 



We have also 



^" 1] (g) = /^ ]t 4 sl (*), (n>*>l). (2-31) 



In terms of the determinant (Wronskian) 



W [A,..., /n](x) d = f det(^M) , (2-32) 

V ClXJ 1 / l<j,k<n 
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(for n = 0, we set W [-](x) = 1.), we have 

W [0 O , 0i, ••• , 05-1, M*) = 0o(x)0i 1] (*)••• ^KWifkaO, (" > s > «)• (2-33) 
Therefore we arrive at the concise formulas due to Crum: 1 ) 
W [0o, 0i, ... , 0«-i](x) = 0o(x)0[ 1] (x) • • • 0fc 1 1] (*) = ± e wW+wW W+ "4^l Wi 



£/W ( x ) = U{x) - 2-^ (log W [0o, 0i, ... , s -i] (*)) , (2-35) 



(2-34) 

cte 2 
, . . . 

W[0 O ,01,... ,(p s -i](x) 
§3. 'Discrete' Quantum Mechanics (pure imaginary shifts) 



4 S] (*) = ^%%11%141M , (n > s > ). (2-36) 



In discrete quantum mechanics, the dynamical variables are, as in ordinary QM, 
the coordinate x, which takes value in an infinite or a semi- infinite or a finite range 
of the real axis and the canonical momentum p, which is realised as a differential 
operator p = —id x . Since the momentum operator appears in exponentiated forms 
e±7P , 7 G 1, in a Hamiltonian, it causes finite pure imaginary shifts in the wavefunc- 
tion e ±7P V( x ) = ip(x=fi'y). This requires the wavefunction as well as other functions 
appearing in the Hamiltonian to be analytic in x within a certain domain including 
the physical region of the coordinate. Let us introduce the *-operation on an analytic 
function, *:/(-►/*. If f(x) = a n x n , a n £ C, then f*(x) d = a* n x n , in which a* 

n n 

is the complex conjugation of a n . Obviously f**(x) = f(x) and f(x)* = f*(x*). If a 
function satisfies /* = /, then it takes real values on the real line. 

The starting point is again a generic one dimensional discrete quantum mechan- 
ics Hamiltonian with discrete semi- infinite energy levels only (|2-ip . Again we assume 
that the groundstate energy is chosen to be zero <?o = 0, so that the Hamiltonian is 
positive semi-definite. The generic factorised Hamiltonian reads 

U = A^A = v / V^e 7P 0^^ + 07^e- 7P v / ^M-^( a; ) ~V*(x), (3-1) 
A = i(ei p ^V^-e~i p ^V(x)), A^ = -i{^/V{x) ei p - y/V*(x)e-&). (3-2) 

By specifying the function V(x), various explicit examples are obtained. 7 ^' 9 )' 11 ) For 
instance, 



%>)(! qi -i,,r 0<*<tt, 0<q<l, 
\a,j\ < 1, {a*, ay, CI3, 04} = {ai, Q2, ^3, 0,4} (as a set), (3-3) 



gives an exactly solvable dynamics whose eigenfunctions consist of the Askey- Wilson 
polynomials 14 ) -16 ) times the groundstate wavefunction 0o(x) (|3-5p . which gives the 
orthogonality measure 0o(x) 2 . 



Crum 's Theorem for 'Discrete ' Quantum Mechanics 



7 



Let us emphasise that the corresponding Schrodinger equation 7iip{x) = £ip(x) 
is a difference equation 

y/V(x)V*(x — ij) ip(x — ij) + \/V*(x)V(x + ij) ip(x + ij) 

- (V(x) + V*(x))ip(x) = Sip(x), (3-4) 

instead of differential in ordinary QM. Again the groundstate wavefunction <f>o(x) is 
determined as a zero mode of A: 

Afa{x) = => H</>o(x) = 0- (3-5) 
The above equation for (fro reads 

^v*{ x - qifaix - q) - ^v{x + q)<h(x + q) = o. (3-e) 

This dictates how the 'phase' of the potential function V is related to that of the 
groundstate wavefunction 4>q. Here we also assume that the groundstate wavefunc- 
tion (j)o(x) has no node and chosen to be real and positive for real x. 

Due to the lack of generic theorems in the theory of difference equations, let us 
assume that all the energy levels are non-degenerate and all the eigenfunctions are 
square-integrable and orthogonal with each other and form a complete basis of the 
Hilbert space: 

H(p n (x) = £ n (p n (x), n = 0,1,2,..., (3-7) 
J 4> n (x)* (j> m {x)dx = h n 5 nm , < h n < oo, n,m = 0,1,2, (3-8) 

In most explicit examples these statements can be verified straightforwardly. For 
simplicity we choose all the eigenfunctions to be real on the real axis (f)^ — <f> n . This 
is possible since H maps a 'real' function to a 'real' function /* = / =>■ (7if)* = Ttf. 

Now the procedure to generate the associated Hamiltonian system goes almost 
parallel with the one shown in the preceding section. We define an associated Hamil- 
tonian T^M by simply changing the order of A and A\ T~t^ = f AA\ as in (|2-lip . 
Then as in (12021) . we have A [0] H [0] = H [1] A [0] and A m H [1] = H [0] A m . Here again 
we have indexed the quantities of the original Hamiltonian system by the superscript 
[°1 , as before. The eigenfunctions and the eigenvalues of the associated Hamiltonian 

are given by 

( x ) d ^ { A [0] ^ ] (x), J <j® (x)*^ (x)dx = £ n h n 5 nm , n, m = 1, 2, . . . , (3-9) 

H [l] (t${x) =£ n <f$(x), ti = 1,2,.... (3-10) 

The same argument as before establishes the iso-spectrality and the one to one 
correspondence of the eigenfunctions of and except for the groundstate. In 
explicit examples given in, 9 ) one can show that <f>'(x) = (f>'(x ; A) = ip(x)/<f>o(x ; A + <5) 
is annihilated by A\ A^(j)'(x) = as in the Darboux transformation. Such functions, 
as before, do not belong to the Hilbert space of the associated Hamiltonian . Here 
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A and 8 are parameters and their shift. For example, tp = 1 for the Meixner-Pollaczek 
case, (f(x) = x for the Wilson case and ip(x) = sinx for the Askey- Wilson case (|3-3p . 
We have used the convention to specify the discrete quantum mechanics by the name 
of the polynomials constituting the main part of the eigenfunctions. 

Thus (/>^(x) is the groundstate wavefunction of with the eigenvalue S\. By 
subtracting it from T6- l \ we can again factorise the positive semi-definite Hamilto- 
nian: 

H [1] =A m A [1] +£!, (3-11) 

A m i( e ip^ v m*(x) - e-% p y/vm (x) ) , 

A [i]\ def e i p - yVM*(x) e~i p ), (3-12) 

with the new potential function V^ 1 ' to be determined now. This imposes quadratic 
relations between the function V"' ' and the unknown V"W: 

V [0] {x - q)V [0] *(x - q) = V [1] (x)V [1] *(x - irf), (3-13) 
V^(x + q) + y [0] *(x - q) = V [1] {x) + vW*(z) - £i, (3-14) 

which are discrete counterparts of the Riccati equation (|2 21|) for the pre-potential W. 
One essential problem is that the connection between the groundstate wavefunction 
<f>o(x) and the function V(x) in the Hamiltonian is indirect, in contrast to the one 
0o (x) = e w W, (E3| in ordinary QM. In discrete QM, the potential function V 
is complex (analytic). The information on the 'absolute value' of V"M(x) can be 
extracted from 

(PEED - The 'phase' part of of V^(x) is given by that of </>^(x) 
through the zero mode equation of A^ (j3- 16[) . just as the 'phase' of 4>o is related to 
that of through the the zero mode equation of A^, (|3-6|) . Thus the following 
formula determining the function vM(x) in terms of the previous function x 
and cj)^ (x) is the main result of the present paper: 

V®(x) d ^ J V W(x-q)VW*{x-il) (3-15) 

n (x) 

A ll Ui ] (x) = 0. (3-16) 

There is no ambiguity in the phase of the square root in the above expression (|3- 15j) . 
It is positive for i = a + i|,oeM ) for which the function inside the square root sign 

is positive definite. Here we have to assume that the groundstate wavefunction (j^ (x) 
has no node, so that the function V"M(x) does not develop unwanted singularities 
in the physical region. It is rather straightforward to verify that (|3-13p is actually 
satisfied by the V' 1 ' in (|315j) . In order to verify the linear relation (|3T4p . one has 
to use the eigenvalue equation for cj)f\ By reversing the order of .A' 1 ^ and A^\ the 
second associated Hamiltonian system Ti} 2 ^ can be defined. Again this process can 
go indefinitely. 
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Here we list the definition of the s-th quantities step by step for s > 1, 

HM^A^A^ + Ss-i, (3-17) 

=>~ 1] ^~ 1] (^ 9 (3-18) 

(x) = Jv^{ x -q)v^*(x-q) ^IS^zM ; (3 . 19) 

0s (a?) 



A [s]t drf _i(^M(a;) e lv _ ^ V [*]*( x ) e -3rP). (3-20) 

Let us note that the phase of the square root in (13T9|) has no ambiguity. Then we 
can show the following for n > s > 0, 

H^{x)=£ n ^(x) 1 (3-21) 

4$* = <f$ : 'real' function, (3-22) 

^W^H( x ) = 0, (3-23) 

ftM =X w U [,1 + £ 4 . (3-24) 

We have also 

4[s-i]t 

<f>t 1] (*)= * _r ^fr), (n>s>l). (3-25) 

The discrete counterpart of the determinant formulas of Crum (|2-33|) - (|2-36|) 
requires a deformation of the Wronskian, the Casorati determinant, which has a 
good limiting property: 

W 7 [/ l5 . . . , f n ](x) ^ det(/ fc (x + i^^-i)) , (3-26) 

V / l<j,k<n 

Iim 7 -3»(«-i)w 7 [/i, / 2 , • . . , f n ](x) = W [/i, / 2 , . . . , / n ](x), (3-27) 
(for n = 0, we set W 7 [-](x) = 1.). Then we have, corresponding to (|2 33|) . 

s-1 

W 7 [0o, 01 , ■ ■ ■ , fa-i, 0n] (x) = J] # ] 0* + *^7) • 4 S] (*) , (3-28) 

fc=0 

<t ] (*) = (3-29) 



Corresponding to (|2-36p . we obtain 

^ W 7 l0o,0l ! ---,0s-lj(a:-*i) 



1=0 
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The proof of the above statements is elementary by induction and the necessary 
nontrivial formulas are only 



Vl s ~ l 1(x + 



a 1 



[X 



q) 



t 1] (x-q) 



(n > s > 1), 

(3-31) 



w 7 [/ , /i, . . . , fs-i, f s ](x + q) w 7 [/ , /i, . . . , fs-i, f n ](x + q) 
w 7 [/ , a, . . . , fs-i, fs](x - q) w 7 [/ , h, . . . , fs-u fn](x - q) 

-iWJfoJi,..., fs-i](x) W 7 [/ Q , /!,... , /,_!, / s , /„](*), (n > s > 0). (3-32) 



§4. Summary and Comments 

Since Crum's paper 1 ) is crisp and elegant, the underlying logical structure is not 
easy to fathom for non-experts or physicists. In section two we reproduce his results 
by a simplest logic in the language of quantum mechanics and using the factorisation 
method, 3 ^ or the so-called super symmetric quantum mechanics, 4 ) so that the simi- 
larity and contrast with the corresponding results of the discrete QM would become 
clear. Due to the lack of essential theorems in the theory of difference equations, 
some important properties of the spectra and eigenfunctions of the generic discrete 
QM must be assumed for the derivation of the associated Hamiltonian systems in 
section three. For example, the hermiticity or the self-adjointness of a discrete QM 
Hamiltonian can only be demonstrated after a proper groundstate wavefunction i^o 
is chosen. 9 ^' 12 ) As mentioned repeatedly, these standard properties are well satisfied 
in explicit examples of discrete QM in. 7 )' 9 )' 11 ) However, to the best of our knowl- 
edge, the general structure of the solutions of the main difference equation (|3-4p as 
well as that for the groundstate (|3-5p for generic potential V(x) and V*(x), has not 
yet been investigated in contradistinction to the ordinary QM. In this connection, 
let us mention two interesting examples of quasi-exactly solvable discrete quantum 
mechanics (example in §IIB1 of the first paper of 12 ) and in §3 of 13 )) in which the 
oscillation theorem does not hold. 

Here are some comments on closely related topics; shape invariance, orthogonal 
polynomials and limiting properties to the ordinary quantum mechanics, etc. 

Shape invariance: The Hamiltonian may contain several parameters A = (Ai, A2, 
. . .) and we write them explicitly TL = H(X), V(x) = V(x ; A), A = -4(A), etc. Let 
us consider the case that the potential function of the first associated Hamiltonian 
VW(x) = ; A) has the same form as the original function V with a different 

set of parameters and up to a multiplicative positive constant k E M + : 



l/W(x;A) = kV(x;X') 



(4-1) 



Here the new set of parameters A' is uniquely determined by A (let us write A' 
si(A)). Then this system has the shape invariance, 5 )' 7 )' 9 ) 



A{\)A(X)^ = /U(A')U(A') +£i(A). 



(4-2) 
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The shape invariance is a sufficient condition for exact solvability. The entire energy 
spectrum and the excited wavefunctions are expressed in terms of £i(A) and 4>o(x ; A) 
as follows: 



n-l 



£ n (X) = £V£ 1 (AM) ) (4-3) 

s=0 

0„(x;A) ocA(X [0] )U(\ [1] )U(X [2] )^--A(X [n ^Mx;X [n] ), (4-4) 
where X [n] is A™ = A, X [n] = si(A [n " 1] ) (n = 1,2, . . .). 

Orthogonal polynomial: Here we consider a generic Hamiltonian (|3-ip of the 

discrete QM. That is the shape invariance is not assumed. Let us define a 'real' 
function r](x) (rf = rj) as a ratio of 4>i(x) and 4>o(x), 

^f\ = a + br,(x), (4-5) 

where a and b (b ^ 0) are real constants. Although rj(x) is not well defined without 
specifying a and 6, this ambiguity (affine transformation of rj(x)) does not affect the 
following discussion. Then (|3- 15[) implies 

Ki'i(. +i j)=n*) fr i7 >7 (1 ' (M) 

ry(xj — 7]{x + 27) 

Let us assume further that the n-th eigenfunction 4>n/4>o is a degree n polynomial in 
this r)(x) for all n > 2: 

n 

^ n (x) = o (x)P n (7/(x)), P n {y) =S ^a rhk y k , V eR > a n ,n / 0. (4-7) 

fc=0 

Obviously ao,o = 1, a = a\$ and 6 = o^i. The orthogonality of the eigenfunctions 
{4> n } implies that {P n (r}(x))} are orthogonal polynomials in rj(x) with respect to the 
weight function 4>o(x) 2 . Then the ratio of the first excited state and the groundstate 
of the s-th associated Hamiltonian system 7iM takes the same form as that of the 
first associated Hamiltonian system TC^ (|4-5|) . 

%M = ^±i£ + ^±W r? M (x ) 5 ^](*) d ^^ + ^7), (4-8) 

4> S (X) °«>* Qs ' S k = Q 

and is related to l/[ s_1 l as in (|4-6p and this goes further down to V^: 

V [s] (x + if 7) = V [s ~ 1] (x + i*=± 7 ) 



4 s ~ l \x + i^ 7 ) - ^-^(x + t*fi 7 ) 
' r?^" 1 ] (x + i^yj - r?! 5 - 1 ] (x + i^ 7 ) 



2 11 1 \^ 1 2 
8 - \W( a; + ^ 7 )- 7/ W( a + i| 7 ) 

V(x)fi f-^-^ + ^\ . (4-9) 
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In all the known examples of exactly solvable QM (ordinary and discrete), in which 
the eigenfunctions have the polynomial form (|4-7|) . the function rj(x) plays a special 
role to define together with the Hamiltonian 7i an algebraic sufficient condition for 
exact solvability. In that case the function t](x) is called the sinusoidal coordinate 
and the sufficient condition is named the closure relation. 8 )> 9 )' n ) As shown above, 
the shape invariance and the closure relation are very closely related. 

Limit from dQM to QM: Here we show that the ordinary QM is obtained 
from the discrete QM in a certain limit by rescaling parameters. Let us introduce a 
positive parameter c and rescale 7 as 7/c and the parameters in V(x) appropriately, 



H = A^A= y/V(x) ecPy/V*(x) +y/V*(x) e~ p y/V{x) - V{x) - V*(x), (4-10) 

A = i(e^y/V*{x) - e~^ p y/V{x)), A 1 = -i[y/V(x) e& - y/V*(x)e-& p ). 

(4-11) 

Assume that V(x) has the following expansion for large c, 

V{x) = a(l + ^w x {x) + 0{^j), (4-12) 

where a is a positive constant. Then we have for large c: 

«£W +0 (i), «tt = * _«*W +0( i), (4.13) 
dx c y&l dx dx c 

d2VVW + 0(I), (444) 





d 


r A 




V«7 


dx 






c 2 

— U -- 


= p 2 + 


aj z 





/ dW(x) 



V dx I dx 2 c 

where the derivative of the pre-potential W(x) is defined by ^^fe) = —Hewi(x). 
Therefore ordinary QM is obtained from discrete QM in the c — > 00 limit. 

Discrete QM with real shifts: Crum's theorem for the discrete QM with real 
shifts can also be formulated in a similar manner. In this case the Hamiltonian is a 
real symmetric tri-diagonal (Jacobi) matrix, either of finite or infinite dimensions. 8 ) 
The factorisation of the positive semi-definite Hamiltonian 7i = A^A also holds 
and A consists of the diagonal and super-diagonal elements only and A^ being its 
transpose, which consists of the diagonal and sub-diagonal elements only. 

Connection with integrable systems: Since the Darboux transformation 2 ) is 
closely related to the inverse scattering method for soliton equations, it is natural to 
ask if the present formulation of Crum's theorem is related to (discrete) integrable 
systems. At present, all the explicit examples considered in discrete QM 7 '' 9 )' 11 ) 
(pure imaginary shifts) have an infinite number of discrete eigenvalues only. In other 
words, the corresponding potentials are confining, that is, they grow to infinity at 
the boundaries or the spatial infinities. Thus the free incoming/outgoing waves at 
infinity do not exist and the corresponding scattering problem cannot be formulated. 
To sum up, we have nothing to report on possible applications of the present Crum's 
theorem to (discrete) integrable systems. 
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After completing this work, we received a recent work by Gaillard and Matveev 23 ^ 
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